Hole Transport in p-Type ZnO 
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A two-band model involving the A- and B -valence bands was adopted to analyze the temperature dependent 
Hall effect measured on N-doped p-type ZnO. The hole transport characteristics (mobilities, and effective Hall 
factor) are calculated using the "relaxation time approximation" as a function of temperature. It is shown that 
the lattice scattering by the acoustic deformation potential is dominant. In the calculation of the scattering rate 
for ionized impurity mechanism, the activation energy of 100 or 170 meV is used at different compensation 
ratios between donor and acceptor concentrations. The theoretical Hall mobility at acceptor concentration of 
7 x 10 18 cm 3 is about 70 cm 2 V s _1 with the activation energy of 100 meV and the compensation ratio of 0.8 
at 300 K. We also found that the compensation ratios conspicuously affected the Hall mobilities. 

PACS numbers: 68.55-a, 71.35c, 78.66.Hf, 81.15.Gh 



I. INTRODUCTION 

Recently, ZnO has received great attention due to its poten- 
tial applications in ultraviolet (UV) and blue optoelectronic 
devices 1,2 and its unique material properties such as a large 
direct band gap of ~3.3 eV at room temperature, high trans- 
parency for visible light, mixability with Mg and Cd to form 
the quaternary system (Zn,Mg,Cd)0 with a tunable band gap 
between 3.0 and 4.0 eV— . Nominally ZnO is grown under 
fairly residual «-type. It is not difficult to achieve higher elec- 
tron concentrations with excess zinc or by doping with Al, Ga, 
orlni^&i. 

On the other hand, it is necessary to obtain p-type layers 
for light-emitter application and make a progress for its reli- 
able technology. But, this has not been solved up to now and 
remains the major bottleneck for ZnO-based optoelectronics. 

Because of the difficulties in growing p-type ZnO films 
as earlier-mentioned, very little is known about the material 
and physical properties of this material. Fortunately, the sit- 
uation is now getting improved and several groups have re- 
ported the growth of p-type ZnO with valid and reasonable 
sets of experimental data. Among them, Look and cowork- 



The calculated hole transport characteristics, with changing 
compensation ratio over a wide range of temperatures (T = 
50 to 400 K), assume transport in both A-valence-hole and 
B-valence-hole bands (a two-band model) with scattering al- 
lowed between these bands. 



II. CALCULATION PROCEDURES 

The analysis follows essentially the same method used for 
representative cubic semiconductors (GaAs, c-GaN), but the 
calculation had to be modified to fit the valence band struc- 
tures of wurtzitic ZnO. In addition, we take the effects of 
warping of surfaces into account for the hole transport stud- 
ied here. The Hall coefficient factor, which generally links the 
experimentally measured Hall mobility with the drift mobility 
can be written as follows: 



< n > 2 



(1) 



ers have reported on the temperature dependence of the Hall- 
k>H | effect measurement. Tsukazaki and coworkers have fabricated 
first ZnO-based light emitters 2 in which they also reported a 
temperature-dependent Hall-effect (T-Hall). In this work, we 
report mainly the theoretical aspects of the T-Hall results and 
the comparison with the experimental data. 

Several groups&ISiAliA^ have already reported the electron 
transport characteristics of «-ZnO both experimentally and 
theoretically. On one hand, there is few theoretical investi- 
gation 1 of the hole transport characteristics of p-ZnO. In this 
paper, the Hall coefficient anisotropy factors (ta) for ZnO are 
determined and theoretical hole transport characteristics (Hall 
and drift mobilities, and effective Hall factor) are calculated. 



where r is the carrier scattering time, < > indicates a ther- 
mal average over the distribution of carrier energy and the in- 
dex i runs over A- (i = 1) and B-valence (i — 2) holes. 

The Hall coefficient factor depends on the degree of warp- 
ing at the warped bands as well as the scattering mechanism, 
and can be written as: 



TAi- 



<rj> 
<n> 2 



(2) 



where rAi is the anisotropy factor for z'-th hole band. 

For their evaluation, one should determine the valence band 
parameters L, M and N defined in terms of the interband ma- 
trix element. Then A, B and C parameters (which are differ- 
ent from the notation related to the splitted valence bands of 
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A 5 


A 6 -A 


—B 


-6.680 


-0.454 6.128 


-2.703 


-2.767 


-4.626 1.41 


1.95 0.796 



TABLE I: Valence band parameters for ZnO. The notation of param- 
eters used here are the same those in ref. 14. 



wurtzitic semiconductors), which have been used in an ap- 
proximate expression for the B-valence- and A-valence-hole 
energy surfaces. We use here the cubic approximation. These 
valence band parameters were calculated from those of effec- 
tive mass parameters (Ai to A e ) recently reported by Fan and 
coworkers^, which are compiled in Table HI The cubic ap- 
proximation yielded the linkage of Ai to A§ with the param- 
eters of L, M and N, whose relations are given by: 



Li =A 2 + A 4 + A 5 , 



L 2 = A u 



Mi = A 2 + A 4 - A 5 



M 2 = Ai+A 3 , 



(3) 



(4) 



(5) 



(6) 



The dispersion relation of energy band is represented for 
the nonparabolic band picture in the following form: 



h 2 k 2 

2m* 



E(1+£E), 



(13) 



where all the notations used take their conventional mean- 
ings. The detailed explanation for the description can be found 
elsewhere^ 

The approximated expressions for the "overlap function" 
describing the hole scattering from the initial state k in band 
i(ki) to the final state k' in band f{k'f) are given by: 



G(ki,k' f ) 



(H 
3(1 



3 cos 2 6) /4 for i = f 
- cos 2 9) /4 for i^f 



intraband 
interband 



(14) 



where 9 is the angle between the wavenumbers of these rele- 
vant two states. 

We summarize the relevant equations related to the scatter- 
ing mechanisms treated in this study. The screening effects for 
the phonon scattering mechanism can be neglected as far as 
the regime of low-impurity concentrations is concerned here. 

The ionized impurity scattering mechanism was treated us- 
ing the Brooks-Herring approximation, where particles inter- 
act via the screened Coulombic potential 1617 . The inverse of 
screening length, A^, can be written for a low hole concentra- 
tion, as follows: 



and 



M 3 = A 2 



N 1 = 2x A 5 , 



N 2 = V2x A 6 , 
where the notations introduced above are: 

L=(L 1 + L 2 )/2, 

M = {Mi + M 2 + M 3 )/3, 
N = (Ni + N 2 )/2 7 



(7) — 



(8) 



(9) 



(10) 



(11) 



(12) 



The relationship between L, M, and N and the well-known 
notations of A, B and C has been given in ref. 1 141 

Here, we used the formulae expressed in ref. 15 to calculate 
Tau the detailed descriptions of which were given in eq. (B.2) 
of ref. Q 

The scattering mechanisms due to (1) deformation potential 
due to acoustic (dp) and optical phonons (npop), polar opti- 
cal phonon (pop), piezoelectric (pe), ionized impurity (imp), 
and interband deformation potentials are taken into account 
in the two-band analysis involving both kinds of hole scatter- 
ings. All scattering rates have been derived assuming a a non- 
parabolic band, which is given below, spherical bands, and 
band warping is accounted for. 



e^knT 



{p+( P + N D )}[l-(N D +p)/N A }} (15) 



where p is the total hole concentration, and e s is the static 
dielectric constant. Other notations have their usual meanings. 
The free hole concentration p can be calculated using equation 
(3.3)ofref.El 

The Hall-effect measurement evaluated the value of E a to 
be approximately 100 meV, while the spectroscopic studies 
experimentally yielded E a ~ 170 meV 8 *^,. Therefore, E a 
used in this work is 100 or 170 meV. 

The scattering rate due to the density of ionized impurities 
Nj (N] = p+2Nu) and the level of ionization of the impurity 
atom Zj is given by: 



Sion(ki, kj) — 



7 



1 



(l-y)G(y) 



Sneooh 3 kf J_i (l-y + l/2fc 2 A 2 ,) 2 



dy, 



(16) 

where y = cos 9 and G(y) — G(cos 9) is the overlap function 
in eq. dl4> . Also, is the effective mass of the relevant band 
into which the holes are being scattered, fcj is the magnitude of 
the initial-state wave vector and the other notations take their 
conventional meanings. 

The scattering rate for the nonpolar optical scattering has 
been given elsewhere 14 . Quantity of D npo is found to be about 
3.9 x 10 11 eV/m. We list here the scattering rate for polar opti- 
cal mechanism because of some possible typographical error 
appearing in ref. [141 
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HI. RESULTS AND DISCUSSION 



is the high-frequency dielectric constant. The 
wavevectors k± correspond to energies Ei ± Ticuq. The equa- 
tions in regard with the scatterings by acoustic deformation 
potential and by the piezo- scattering have been enumerated 
in ref. 14 as equations (c.2) and (c.3). 

In order to determine drift and Hall mobilities from the 
above-mentioned scattering mechanisms, we have determined 
the effective Hall factor using the "relaxation time approxima- 
tion" or otherwise known as "Mathiessen's rule". The effec- 
tive Hall factor r e g is: 



r s 



(1 + a 3 / 2 )(r ia 3 / 2 + (3 2 r 2 ) 
(l + a 3 / 2 /3) 2 



(18) 



with, a — mi — rri2, [3 = /ii///2> Mi = e < Tj > /m*, 
and Ti = tai < rf > / < Ti > 2 . Quantities m, /Xj, and 
< tj > are the effective masses, drift mobilities and the re- 
laxation times for the z'-th hole band, respectively. (The index 
i = 1 always refers to A-valence, and the index i = 2 to 
B-valence holes.) 

The drift mobility p, is determined from relationships be- 
tween drift and Hall mobilities given by: 



r c s' 



(19) 



where fijj is the Hall mobility and r e g is the effective Hall 
factor given by eq. dl 8I >. 

The Hall mobility, /.in, is determined to be: 



Ph 



e r A \ <t{> 



< n > +U- 1 /' 2 < T 2 > 



(20) 



where to* is the effective A-valence-hole mass, and tax and 
t A2 are given in eq. (B.2) of ref.ll4l 



Property 



Units Value 



Eg Energy gap for A-valence hole 

E g Energy gap for B-valence hole 

mi /mo Effective A-valence hole mass ratio 

w-2/mo Effective B-valence hole mass ratio 

p Density 

s Sound velocity 

Ei Acoustic deformation potential 

Dnpo Optical coupling constant 

e s Static dielectric constant 

£00 High-frequency dielectric constant 

huio Polar optical phonon energy 

h„z Piezoelectric constant 



cV 
cV 



kg/iri 
m/s 
eV 

eV/m 



eV 
C/m 2 



3.445 
3.455 
0.59 
0.59 
5.67526 
6.09 x10 s 

15.0 
3.9xlO n 
8.2 
3.7 
0.072 
0.89 



TABLE II: Material parameters of ZnO used in the calculation. 

The material parameters used in the calculations on ZnO are 
listed in TablelTTl Rode 9 suggests a theoretical value of 3.8 eV 
for deformation potential Ex, but the later studies resulted in a 
better agreement with the experimental data if the larger value 
was used. 



We determined the valence band parameters (L, M, N, A, 
B, C) using Luttinger parameters. Their results are listed in 
Table I. The evaluated anisotropy factors tax and ta2 both 
yielded the values very close to the unity, which is signifi- 
cantly different from the situation observed in GaN. This dif- 
ference may be due to the isotropic nature of the hole masses 
in ZnQA2. 
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FIG. 1: Variation with temperature of calculated scattering relax- 
ation time corresponding to nonpolar optical phonon mode (a), po- 
lar optical phonon mode (b), acoustic deformation potential mode 
(c) and piezoelectric mode (d). Many traces are almost overlapped 
each other because the two hole masses are the same in ZnO. 
(h^h - scattering from A-valence-to-A-valence-hole band, h^l - 
scattering from A-valence-to-B-valence-hole band, 1— >1 - scattering 
from B-valence-to-B-valence-hole band, and 1— »h- scattering from 
B-valence-to- A- valence-hole band.) 

The partial '/u versus temperature' curves are plotted in 
Fig-E The interband scattering characteristic due to nonpolar 
optical phonons shown in Fig. ^a) is different from that due 
to other scattering mechanisms shown in Figs. nW,Ec) and 
^d). This is because the scattering relaxation time due to non- 
polar optical phonon scattering is determined by the final-state 
effective hole mass (m%), while that due to other scattering 
mechanisms is determined by both m*f and the magnitude of 
the initial state of the wave vector (fej) in which nonparabol- 
icity is considered. We remark that, the relative importance 
of the various scattering mechanisms mimics that found in 
ZnSe 20 . For example, at temperatures below 300 K, the acous- 
tic deformation potential scattering plays the most important 
role. 
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FIG. 2: Variation of calculated average scattering time for ionized 
impurity scattering with temperature for E a = 100 and 170 meV 
with dependence on compensation ratio {K/Nd = Na)- The accep- 
tor concentration (Na) was set to 7 x 10 18 cm -3 in the calculations. 



Fig.[2]plotted average scattering time (r av ) by ionized im- 
purities calculated for two activation energies (E a = 100 and 
170 meV) with compensation ratios (K = Np /Na) between 
0.1 and 0.8. Value of r av is calculated with r s ~^ t as the sum 
of the scattering rate from A-to A-hole, A-to-B-hole, B-to-B- 
hole and B to- A-hole bands. 

In the range of those acceptor activation energies (E a — 
100 — 170 meV), the effect of E a on the average 
scattering time at high temperatures is fairly negligible 
(r ol) (170)/r ot ,(100) ~ 1 at 300 K), while the corresponding 
effect for a low compensation ratio (K = 0.1) is appreciable 
(r a „(170)/r a „(100)~2at 300 K). 
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FIG. 3: Variation of calculated effective Hall factors with tempera- 
ture due to lattice scattering mechanisms for A- valence holes (a) and 
for B-valence holes (b). npop is the nonpolar optical phonon scat- 
tering, pop is the polar optical phonon scattering, dp is the acoustic 
phonon scattering, pe is the piezoelectric scattering, total is the total 
of all four lattice scattering mechanisms. 

The temperature dependences of Hall factors for the four 
different scattering mechanisms plotted in Fig. [3] The Hall 
factors for A-valence-hole and B-valence-hole bands of non- 
polar and polar optical scattering vary with temperature, while 



those of the acoustic deformation potential and piezoelectric 
scattering remained nearly unchanged, the different observa- 
tions of which may be attributed to the temperature depen- 
dence of the scattering relaxation time of respective partici- 
pating phonons. The values of the Hall factor both for A- 
valence-hole bands and for B-valence-hole bands varied from 
1.0 to 1.5. 
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FIG. 4: Theoretical drift mobilities due to each lattice scattering 
mechanisms (a), and the theoretical Hall and drift mobilities and the 
effective Hall factor (b) due to lattice scattering mechanisms as a 
function of temperature. 

The partial Hall and drift mobilities, and the effective Hall 
factor that we calculated are plotted in Fig.|4]against temper- 
ature. The acoustic phonon scattering is the most important 
mechanism, limiting the hole mobility over a wide range of 
temperature. For temperatures above 300 K, the polar optical 
phonon scattering mechanism is the dominant factor. 

The obtained effective Hall factors vary from 1 . 1 to 1 .7 on 
changing the temperature, as shown in FigJUb), indicating the 
importance of proper considerations of the Hall factors at high 
temperatures when comparing the drift mobility with fijj- 

Finally, Fig. |5] corresponds to the similar plot with Fig.0] 
including now the ionized impurity scattering with four differ- 
ent compensation ratios (K) between 0.0 and 0.8. The accep- 
tor concentration is Na = 7 x 10 17 cm -3 and E a = 100 meV. 
We also plotted the experimental data of p-ZnO with nitrogen 
concentration of 2 x 10 20 cm -3 (different from the accep- 
tor concentration Na) for comparison. The compensation ra- 
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FIG. 5: Theoretical Hall mobilities (a), and effective Hall factor (b) 
due to lattice and ionized impurity scattering mechanisms for differ- 
ent ratio of compensation as a function of temperature. E a used in 
the ionized impurity scattering mechanism is 100 meV. Also shown 
in (a) is a set of experimental data of ZnO:N thin film for comparison. 



tio and free hole concentration at 330 K of this sample are 
p w 8 x 10 16 cm -3 and approximately K = 0.8^. 

Unfortunately, the calculated Hall mobility is still higher 
than that in the experimental reports. We considered only dp, 
npop, pop, and imp scattering mechanisms in our modeling, 
which may be a result of this poor agreement. A dislocation 
scattering and the space-charge effect probably also contribute 
as a scattering mechanism of holes to the mobility determina- 
tion. Actually, it is known that the former mechanism have 
an influence on relatively-poor-quality «-type ZnO epitaxial 
layers^. We would like to propose more experiments using 
further optimized samples be done in future. 

With an increaase in the compensation ratio {K), the peak- 
ing temperature at maximum mobility also increases. The ef- 
fective Hall factor (r e ff) increases as K decreases with increas- 
ing temperature. From Figs.|5Ja) and|5Jb), it is observed that 



the maximum value of the Hall and drift mobilities are in the 
temperature range 50 to 100 K for all values of K between 0. 1 
and 0.8. In comparison with Fig. 0Jb), Fig.|5jb) shows that the 
inclusion of scattering by ionized impurities into the Hall fac- 
tor calculation is important at high compensation ratios (K = 
0.1 to 0.8), exhibiting wavy behavior with the change in tem- 
perature. The obtained effective Hall factors vary from 1 . 1 to 
1.5 with changing temperature (K > 0.1). The temperature 
at which the effective Hall factor is minimized and the effec- 
tive Hall factor is affected by the ionized impurity scattering 
mechanism. 



IV. CONCLUSIONS 

In this study, a two-band model analysis was carried out in 
detail on nitrogen-doped ZnO epitaxial layers. The theoretical 
Hall and drift mobilities, and the effective Hall factor were 
calculated at temperatures of 50 to 400 K. The hole transport 
is characterized by the scattering by the acoustic deformation 
potentials at temperatures below 330 K. The theoretical Hall 
mobility at acceptor concentration Na = 7 x 10 18 cm~ 3 is 
about 70 cm 2 /V _1 s —1 with the activation energy of 100 meV 
and the compensation ratio of 0.8 at 300 K. The compensation 
ratio strongly affected the hole transport characteristics. 
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